STUDYING THE TOPOLOGICAL STABILITY OF THE A$ 4 KINK. 
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The lambda-phi4 kink is linearly and topologically stable. We study how extra energy perturbations are 
dissipated beyond the linear regime. We found that depending on the width, amplitude and energy of a Gaussian 
perturbation different scenarios are possible: radiation, oscillons creation, kink anti-kink pairs production and 
shock waves. 
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1 Basic information 

Scalar fields with a A$ 4 potentials has been studied along several decades, it has been used to model the behavior 
of physical systems in several branches of physics like condensed matter, plasma physics and field theory among 
others. This toy model is interesting because although it is one of the of simplest models, it is non-integrable [T|. 
This model has a kink and a anti-kink solution, see pj]. The kink is linear and topologically stable, meaning that it 
is a kind of 'indestructible' solution of the field equations. In the mid 70's it was found numerically, that after the 
collision of a kink with an antikink, something like a breather solution appeared, see [3]. Segur and Kruskal showed 
that a breather is unable to exist [3] , but it is an unstable oscillatory solution of the of the A<1> 4 theory, that has 
an unusual long time existence. After that article, new work have been done studying this oscillon, its existence, 
stability (see for instance [5]) and its interaction with kinks. 

The aim of the present paper is to revisit the behavior of the of the A$ 4 kink under perturbations. We know 
that it is stable, but we want to know how extra energy of a perturbation is dissipated. It is well known that when 
the perturbation is small enough, the extra energy is dissipated through radiation. Our interest is to study beyond 
this regime, increasing the amplitude beyond the linear domain and further if possible. 
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2 The model 




(i) 



(2) 



We use units such that c = 1. This equation has the static solutions, 




(3) 



that are known as kinks (+) and anti-kinks (— ). 
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We shall perturb the kink with a perturbation of Gaussian profile. This profile is chosen because it is smooth and 
localized. In this way, our perturbation problem can be defined as solving equation ^ with the initial condition, 

0(x,O) = tanh (^x) + Ae~ b{ - x - Xo)2 , (4) 

where the constants A, b and Xq describe the amplitude, inverse of width and position of the Gaussian perturbation. 



3 The method and preliminary results 

To solve our differential equation we use a finite difference scheme, for simplicity we use the Leap-Frog scheme that 
is of the order two in time and space and it is also explicit, then it is not difficult to implement. Fixed boundaries 
far enough from the kink and the perturbation are imposed. We note that three parameters can be tuned in the 
perturbation A, b and Xq. 

Different results are obtained depending on the amplitude (A) and the width (1/6). When b is of the order of 
0.1 and the amplitude is less than one, basically all the perturbations are transformed in radiation. This behavior 
is shown in figure [l] When we further increase the amplitude, an oscillon is created, see figure [2j Beyond the 
threshold of four times the energy of the kink, it can appear a kink-anti kink pair and if the energy is high enough 
an oscillon, see figure [3] 

If the amplitude is further increased, the double well is ignored and only the A<I> 4 potential is seen. In this case 
it is observed oscillations around = 0. See for example figure [4j where A = 3 and the total enery of the system is 
almost 70 times the energy of the kink. Then, the perturbation is spread, its amplitude decreases while going away 
from the kink. When the field is of the same order of unity, creation of structures occurs, see figure [5] 

If b is smaller, b ~ 0.01, i.e. the perturbation is broader similar results are obtained. But the formation of 
oscillons at lower energies is more likely. For b high, b ~ 100, the perturbation is very sharp and it seems that the 
energy goes away from the kink like shock waves, see figure [6j 

It is very interesting the different behaviors of the field under different perturbations. In figure [7] we show the 
parameter space of the perturbation. The curves are equal energy lines, the legend shows the value in terms of 
the kink energy, i.e. 2 means total energy twice the kink. Also this figure shows the outcome of the perturbation 
at the time t = 100 units, "rad.", means radiation, "osc." an oscillon was created, "k+ak" a kink anti-kink pair, 
'"k+ak+st" a kink anti-kink plus an oscillon or two, "phi" means that the double well is unnoticed and finally 
"shock" stands for a lot energy is in the border, like figure [6] For a given amplitud we can obtain the same energy 
with two widths. If is sharp enough, the spatial derivative term in the kinetic energy is huge, and on the opposite 
case we have a lot of energy because of the bigger area of the perturbation. 

4 Discussion 

We did not find in the literature where large perturbations were explored. Due to nonlinear character of the field 
equation, some precautions have to be taken, some regions will oscillate rather quickly in time and space and it is 
important to have a small grid spacing. But, there are zones where this small grid spacing is not needed. There 
are schema like the AMR (Adaptive Mesh Refinement), where the refinement is applied intelligently to adjust itself 
and reduce calculations. 
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Perturbation with b=0.1, A=0.5 and time t=100 
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Figure 1: With a small initial amplitude, all energy dissipates in radiation. 



Perturbation with b=0.1. A=l.l and time t=100 
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Figure 2: If the amplitude is increased, an oscillon is formed. 
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Perturbation with b=0.1, A=1.5 and time t=100 
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Figure 3: When energies go over four times the kink energy, kink-anti kinks can be created. 



Perturbation with b=0.1, A=3 and time t= 10 
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Figure 4: Almost 70 times the kink enery, it shows oscillation around <j> = 0, the double well is unnoticed. 
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Perturbation with b=0.1, A=3 and time t=200 
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Figure 5: Almost 70 times the kink enery, Some time later, the creation of structures begin. 
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Perturbation with b= 100, A=2 and time t= 100 
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Figure 6: When the perturbation is very sharp, b ~ 100. 
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Figure 7: When the perturbation is very sharp, b ~ 100. 
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